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Figure  captions. 
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fig  1:  gas  bubble  around  a  crack  tip. 


fig  2:  tension  regions  for  a  gas  bubble  around  a  crack  tip. 


fig  3:  mode-1  crack  with  a  line  distribtuion  of  edge  dislocations. 


fig  4:  mode-1  crack  with  two  symmetric  line  distribtuions  of  edge  dislocations. 


Fig  5:  Plastic  zone  with  distance  of  the  gas  bubble  from  the  crack  tip  (a  =  55°). 


Fig  6:  Plastic  zone  with  distance  of  the  gas  bubble  from  the  crack  tip  (a  =  60"). 


Fig  7:  Plastic  zone  with  distance  of  the  gas  bubble  from  the  crack  tip  (a  =  65°). 


Fig  S:  Dislocation  density  with  and  without  gas  bubble. 


Fig  9:  Plastic  zone  for  two  symmetric  dislocation  lines  with  and  without  gas  bubble. 


Fig  10:  plastic  zone  with  the  number  of  gas  bubbles. 
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Abstract 


A 

During  the  loading  of  a  crack  edge  dislocations  are  emitted  from  the  crack  tip 
forming  a  plastic  zone.  The  presence  of  gas  bubbles  near  the  crack  tip  can  change 
the  plastic  zone  size  and  location  drastically.  A  detailed  study  of  the  force  field 
on  an  edge  dislocation  due  to  the  crack  load  and  gas  bubbles  allows  for  the 
determination  of  the  plastic  zone  as  a  function  of  the  inclination  angle.  The 
dislocation  distribution  along  the  plastic  zone  is  found  from  force  equilibrium 
considerations.  Due  to  the  presence  of  gas  bubbles  a  shrinking  and  repulsion  of 
the  plastic  zone  from  the  crack  tip  is  noticed.  This  process  indicates  the  material 
embrittlement.  The  dislocation  free  zone  (DFZ)  near  the  crack  tip  grows.  This 
embrittlement  can  become  crucial  in  material  design  for  fusion  applications. 
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Introduction 


The  degradation  of  the  inner  walls  of  fusion  reactors  is  a  serious  problem.  The 
bombarding  particles  cause  nuclear  reactions  which  may  produce  considerable 
concentration  of  foreign  elements  inside  the  material  (Ullmaier  1984).  Here  we 
concentrate  on  the  (n,a)  reactions  which  generate  helium  gas  within  the  metals. 
Because  of  the  low  diffusion  rate  the  helium  basically  stays  where  it  is  formed  and 
can  exert  there  a  considerable  pressure  (Haubold  1983  and  Jaeger  et  al.  1983). 

In  this  article  we  study  the  influence  of  this  bubble  formation  on  the  embrittlement 
of  the  material.  Experiments  have  shown  that  plastic  zones  are  formed  during 
crack  propagation  due  to  the  emission  of  dislocations.  Several  investigators  have 
discussed  the  pilc-up  of  screw  dislocations  in  thin  regions  near  the  crack  tip 
forming  a  plastic  zone  (Chang  et  al.  1981  and  1984).  Here  we  develope  a  plasticity 
model  for  the  thicker  areas.  Observations  by  means  of  the  electron  microscope 
have  shown  that  edge  dislocations  are  emitted  occupying  regions  within  angles  45 
to  90  °  to  the  crack  plane  causing  a  zigzag  propagation  of  the  crack  (Horton  et 
al.  19S2).  Our  model  is  able  to  point  out  to  the  location  and  the  size  of  the  plastic 
zone  for  different  angles  of  inclination. 

Employing  the  present  model  we  would  like  to  explain  the  so-called  embrittlement 
in  materials  containing  helium  gas.  Since  all  materials  contain  microcracks,  under 
external  loads  dislocations  are  emitted  from  the  crack  tip  and  plastic  zones  are 
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■  formed.  Our  model  shows  that  the  presence  of  helium  bubbles  reduces  the  plastic 

) 

|  zone  size.  This  we  believe  is  one  of  the  prominent  mechanisms  underlying  the 

process  of  embrittlement. 

I  Stress  Field  in  a  plate  with  a  line  crack,  under  a 

pressure  source. 

In  order  to  determine  the  force  field  on  an  edge  dislocation  due  to  a  gas  bubble 
in  the  neighbourhood  of  a  crack,  we  need  to  know  the  stress  field.  To  solve  this 
problem  wc  use  the  superposition  principle.  Consider  aline  crack  |x|  <  c,y  =  z  -  0, 
located  in  the  plane  z  =  0  of  rectangular  coordinates  (x,y,z)  and  a  gas  bubble  at  a 
point  (l,  t])  outside  the  CTack,  fig  1.  The  pressure  of  the  gas  bubble  will  generate 
stresses  on  the  crack  surfaces.  To  free  the  crack  surfaces  from  tractions  we  super¬ 
pose  the  negative  of  the  tractions  induced  by  the  bubble  on  the  track  surfaces. 
The  effect  of  the  stress  induced  on  the  surface  of  the  gas  bubble  due  to  this 
superposition  is  in  the  neighborhood  of  0. 1  %  of  the  effect  of  the  original  bubble 
pressure  for  the  configuration  studied.  Hence  we  will  neglect  this  effect. 

The  total  stress  field  is  given  by: 

a‘kl(x.y)  =  Oki(x.y)  +  °£(*,y)  k,l  =  x,y  (1) 


Stress  Field  in  a  plate  with  a  line  crack,  under  a  pressure  source. 
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where 

c'ki(x <}')'■  total  stress  field 

g£/(jc,  j):  stress  field  due  to  the  gas  bubble 

c^i(x,y):  stress  field  due  to  the  interaction  of  the  gas  bubble  with  the  crack. 


It  is  important  to  note  that  the  assumption  of  stress-free  crack  surfaces  is  only 
valid  if  the  crack  opens  up,  i.e.  if  the  gas  bubble  generates  tension  along  the  crack. 
For  a  gas  bubble  located  at  the  origin  of  the  coordinates,  the  normal  stress  oyy  is 
given  by  (Erincen  19S0): 


2  *2  -  y2 

°y>  =  aP“ - ^ 

(*  +  y2) 


(2) 


where  a  is  the  bubble  radius  and  p  is  the  bubble  pressure. In  order  to  generate  a 
tensile  load  on  the  crack  surface  the  bubble  must  be  confined  to  a  region  of  ±  45° 
from  the  line  of  the  crack,  at  the  crack  tip  (fig  2). 


We  assume  that  the  crack  is  symmetrically  loaded,  i.e.  for  a  bubble  situated  at 
(b,Ti)  there  is  a  symmetrically  situated  bubble  at  (-£,ii)  left  of  the  y-axis.  This 
results  in  a  symmetric  tensile  stress  and  an  antisymmetric  shear  stress.  For  our 
configurations  the  interaction  between  the  two  bubbles  is  negligible  (Ling  1948). 


In  this  way  we  split  up  the  problem  in  two  easily  solved  sub-units: 
1.  the  stress  due  to  a  gas  bubble  is  given  by: 


2  (x  -  l)2  ~  (y  -  tl)2 

3  '  [(*  -  4)2  +  O'  ~  ti)2]2 


(3) 


Stress  Field  in  a  plate  with  a  line  crack,  under  a  pressure  source. 
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(4) 
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°yy  “  °xx 


°xy  =  "  ^‘P 


2  (*  ~  E)(y  -  ii) 


[(*  -  l)7  +  (y  -  h)2] 


2.  The  calculation  of  the  stress  field  due  to  an  arbitrary  symmetrically  loaded  crack 
reduces  to  a  superposition  of  a  symmetric  mode-1  problem  and  an  antisymmetric 
mode-11  problem. 


i. 

cM 


The  Boundary  Conditions  (B.C.)  of  the  symmetric  mode  I  problem  are: 


°xyly=C  =  0 


1 

Vi 

I 

•i§ 
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rJjj 
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ijy 

‘V? 

ijS 

\yjm 

♦S] 


ayylcrack  “  P(x) 


^•loutsidc  crack  ® 


p(x)  =  p(  -  x) 


where  ckl,  (k,l  =  x,y)  denote  the  stress  tensor,  (ux,uy)  the  displacement  vector 
and  p(x)  is  the  traction  applied  to  the  crack  surface. 


The  stress  field  for  this  problem  is  given  by  (Sneddon  et  al.  1969,  p  26): 


°xx(x,y)  =  v'Ii0  (y*  "  l)*V(S)e  4ycos(£x)d£ 


°yj(*.y)=  (>^  +  !)*v(*)e  ^ycos(£x)d4 
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Stress  Field  in  a  plate  with  a  line  crack,  under  a  pressure  source. 


ox%(x,y)  =  ->.  f  J  £2v(4)c  ~  ^  sin(4x)d^ 

'  '  o 

(12) 

o„(x.  y)  =  Vfo^  +  Oyy) 

(13) 

where 

v(=)  =  N  §  Jo‘  f(t)J0(£t)dt 

(14) 

,  ,  2t  r-  P(u)du 

fl‘>  "  n  Jf  - - t 

‘  V  -  U3 

(15) 

Here  Jn(:)  denotes  the  Bessel's  function  of  the  first  kind  of  order  n. 

The  B.C.  of  the  antisymmetric  mode  11  problem  are: 

O 

II 

o 

II 

£ 
t r 

(16) 

°xylcrack  q(x) 

(17) 

■  outside  craek  —  0 

(18) 

q(x)  =  -q(-x) 

(19) 

where  q(x)  is  the  shear  stress  applied  to  the  crack  surface. 

Stress  Field  in  a  plate  with  a  line  crack,  under  a  pressure  source. 


The  solution  of  this  problem  is: 


=  “  \  n  i  “  2);y(;)e  ”'cos(;x)di 


(20) 


o>7(x,y)  =  v  £  jo  42>v(;)e  ^'cos(4x)d4 
axy(*’>')  =  N'  i  J  (“V  "  l);v(;)e“^  sin(£x)d' 


(21) 


(22) 


^2z(^’  ^  +  ^yy) 


(23) 
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where 


o^(x.y)  =  J  U0(n-t)[yh(x,y,t)  -  f(x,y,t)]dt 


°vy(x- >)  =  “J  Uc(nrt)[ux,y,t) +yh(x,y,t)]dt 


o"(x,y)=  -yj  tJ0(nnt)j(x, y, t)dt 


Antisymmetric  mode-II: 


q(x)  =  2_  qnsin(n"*)  +  qcx 


°k:  =  L  qr.ak!(x'  >')  (k,  1  =  x,  y) 


where 


a^*(x,y)  =  -  J'  g.(t)[yh*(x,y,t)  -  2f* (x,y,t)]dt 


o"  (x,y)  =  yj  gn(t)h*(x,y,t)dt 


ojv(x.y)  =  J1gn(t)[yhs(x,  y,t)  -  fs(x,y,t)]dt 


gn(t)  =  tJ](nnt)  n  ^  1 


goto  =  t2  /2 


The  term  q0x  in  (31)  is  to  account  for  a  non  zero  shear  stress  at  the  crack  tip. 
The  functions  J[x,y,t),  h{x,y,t).  j{x,y,i),  f  (x,y,t),  h'{x,y,  t),  fs{x,y,t)  and 
Aj(x,y,f)  are  collected  in  the  appendix. 


Stress  Field  in  a  plate  with  a  line  crack,  under  a  pressure  source. 


Along  the  crack  line  a  great  deal  of  simplifications  are  achieved: 


Symmetric  mode-1 


f,  U0(n— t)dt 

(x  -  t  ) 

(38) 

0)  =  o^(x,  0) 

(39) 

o^(x.m  =  o 

(40) 

Antisymmetric  mode-11 

o^(x,  0)  =  0 

(41) 

Q 

p 

o 

II 

o 

(42) 

n  r1  tf(t)dt 

Ox).(x,  0)  -  J  32 

■  (X2  -  tV 

(43) 

Method  of  approach 


Experimental  observations  indicate  that  from  the  crack  tip  emanate  one,  or  two 
sets  of  edge  dislocations  making  an  average  angle  ±a  with  the  crack  line.  We 
idealize  this  picture  by  assuming  that  the  edge  dislocations  constitute  two  straight 
lines. 

Z  and  t  are  local  coordinates  along  the  plastic  zone  (t,  ^e[0,Z.]),  fig  1.  £2,  C3  and 
are  symmetrically  located  points  with  respect  to  £:  if  £  corresponds  to  (x,y), 
then  C;.  corresponds  to  (-x,y),  to  (-x,-y)  and  ^  to  (x,-y). 

At  first  we  look  at  the  case  of  just  one  set  of  dislocations  (fig  3)  at  each  crack  tip. 

We  distinguish  five  kinds  of  forces  on  any  of  these  dislocations: 

1.  The  force  F e  due  to  all  'external'  sources,  i.e  sources  not  linked  to  the  dislo¬ 
cations  generated. 

Fe(t)  =  b(t)F0(t)  (46) 

where  b(t)  is  the  Burger's  vector  density  at  t, 


y> 

$ 


•'yw 


F0(l)  =  {[cxy(i)cos(a)  +  OjT(t)sin(a)]'e; 


-  [o^t'l  cos(a)  +  o^(t)  sin(a)]T 


o{;[t)  (k.l  =  x.y)  is  the  'external'  stress  field  at  t. 


In  the  case  under  investieation  we  have: 


°ki(t)  =  Oki(0  +  aid(l)  +  °kl(l) 


and  hence  o£/  consists  of  three  contributions: 


oCki(i):  the  stress  field  due  to  the  external  crack  loading. 


the  stress  field  due  to  the  gas  bubbles. 

chK)(t):  the  stress  field  due  to  the  interaction  of  the  gas  bubbles  with  the  crack. 


2.  The  force  F  d  due  to  the  other  dislocations  in  zone  1 . 


Fd(t)  =  b(t)JLK1(t,;)bQd; 


where 


K,(t.C)  “ 


2n(I  -  v)(t  -  5) 


Here  p  is  the  shear  modulus  and  v  is  Poisson's  ratio. 


Method  of  approach. 


3.  The  force  F &  due  to  the  y-symmetrically  located  dislocations. 


hds(t)  =  b(t)  P  bQK;(t,:)d; 


(51) 


where 


K2(t,;)  =  |[o^.(t,;)cos(a)  +  Oj^(t,  C)  sin(a)] 


-  [°x\(t-  ^)cos(a)  -r  o^(t,;)sin(a)]‘e  y| 


(52) 


o is  the  stress  Geld  at  t  due  to  a  unit  Dirac-delta  dislocation  density  at  £2- 


4.  The  force  F  cj  due  to  the  interaction  of  all  dislocations  with  the  crack. 


FaJ(0  =  b(t)JLbQK3(t,;)d; 


(53) 


where 


K  3(t, ;)  =  {[o^.(t,;)cos(a)  +  og(t,;)sin(a)]' 


-  [a^(t,;)cos(a)  +  o^.(t,;)sin(a)]'e  y | 


(54) 


a *f(f,  £)  Is  'be  stress  field  at  t  due  to  the  interaction  of  the  crack  with  a  unit 


Dirac-delta  dislocation  density  at  £  and  at 


5.  The  force  F  j  due  to  friction  in  the  glide  plane  of  the  dislocation. 


F(<t)  =  -Cfb(t)  c  t 


(55) 


Method  of  approach. 
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where  <y  is  of  the  order  of 


cj-  =  10  "V  a  10  2|i 


(56) 


The  equilibrium  configuration  of  the  system  requires  that: 


Fe  +  rd  +  Ff  =  0 


(57) 


or 


Fo(»)  +  JoL  [K3(t.;)  +  K2(t,;)+  K3(t,;)]b(;)d;-cfTl  =  0  (58) 


At  first  we  neglect  the  interaction  forces  F  &  and  F Cd-  Hence  we  have: 


F0(t)  +  JLK,(t-;)b(;)d;-cf'el=  0 


(59) 


Due  to  the  chosen  configuration  the  forces  between  the  dislocations  F^are  parallel 
to  their  line  of  position,  i.e.  they  have  the  same  angle  of  inclination.  Indeed, 
substituting  (50)  in  (59)  one  obtains: 


F0(t)  +  et[joLK,(t,;)b(;)d;-cf]  =  o 


(60) 


where 


K,(t.;)  =  iK,(t.;)i 


(61) 


Method  of  approach. 
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Hence,  in  order  to  maintain  the  equilibrium  the  forces  due  to  the  crack  load  must 
also  have  an  inclination  angle  a  (a  necessary  condition).  From  this  we  deduce  the 
following  approach  to  the  solution  of  our  problem: 

i.  We  choose  an  angle  a. 

ii.  We  determine  the  force  field  F c  due  to  the  crack  load  on  an  edge  dislocation 
with  this  angle  of  inclination. 

iii.  In  this  force  field  we  look  for  a  straight  line  with  inclination  angle  a,  of  forces 
having  this  same  inclination  angle.  This  gives  us  the  size  and  position  of  the  plastic 
zone. 

iv.  Knowing  the  fo'cc  field  along  this  line,  we  can  calculate  the  edge  dislocation 
distribution  (Burger's  vector)  which  will  generate  a  force  field  with  the  opposite 
sign  (force  equilibrium).  This  requires  the  solution  of  the  following  integral  equation: 

jLK,(t,;)b(;id;= -f(t)  (62) 

o 


where 


f(l)  -  Fg(t)  cf 

(63) 

F  o(*)  =  1  F  oWI 

(64) 

v.  The  interaction  forces  betw  een  crack  and  dislocations  can  be  treated  by  iteration. 
Let  the  solution  of  (62)  be  bo(i).  The  dislocation  density  b\(i),  after  one  iteration 
is  given  by  the  solution  of 


?0(t)  +  £  k  ,(t.;)b,(;)d;  +  JoL  [k2(i,;)  +  K3(t,;)]b0(;)d;  -  Cf‘?t  =  o  (65) 


Method  of  approach. 
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The  iterated  location  of  the  plastic  zone  is  determined  by  the  force  field 


I  |(t)  I  0^1  ""  j  [^jtt.,)  K.  j(t,  ,)  jby(,'d. 


(66) 


and  the  density  b-(i)  is  the  solution  of 


fLK,(t,;ib,(;id;  =  -f,(t) 

where 


(67) 


f,(t)  =  |F,i  -  cf  (68) 

In  case  of  two  symmetrically  located  dislocation  lines  (fig  4),  our  equations  do 
not  change  on  condition  that  we  redefine  o*;  and  o£f: 

')  is  the  stress  field  at  t  due  to  a  unit  Dirac-delta  dislocation  density  at 
^3  and  i.c.  all  x-  an  y-  symmetric  positions  of  t  = 

o^;(r.  is  the  stress  field  at  t  due  to  the  interaction  of  the  crack  with  a  unit 
Dirac-delta  dislocation  density  at  t  =  £  and  all  its  x-  and  y-symmetric  positions. 


Due  to  the  significance  of  the  correction  forces,  iteration  is  more  important  in  the 
symmetric  case. 


The  determination  of  the  edge  dislocation  distribution  corresponding  to  the  cal¬ 
culated  force  field  amounts  to  the  solution  of  the  following  Cauchy  integral  equa¬ 
tion  for  b(x)  (62): 


ft  bQdd; 

«  2n(  1  -  v)(;  -  t) 


-ftt) 


(69) 


Method  of  approach 
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(70) 


1  =  f(.t)  0  £t  £  L 

■o  ,  -  t 


where 


mMC) 


(71) 


The  solution  of  (70)  which  is  zero  at  L  is  given  by  Nlikhlin  1979  p  130: 


— { f  "L  ^  KK-rMrl  fQd;|  (72) 
7t\  t(i.  - 1)  C  1  1  “  ’  »  LL  -  ’J  ) 


1  2 


C*(t)  = 


te[0,L] 


(73) 


The  first  integral  is  a  Cauch>  integral  whose  singularity  can  be  removed  as  follows: 
we  write  the  integral  as  one  centered  about  the  singularity  and  the  rest:  for  i  £  L  2: 


C  t  - 


.1  2 


f  l(t  -  u)(L  -  t  +  u)]'  2f(t  -  u)^  +  JL1^— yl  fU)d; 

_  l  2t  I 


(74) 


Then  w  e  subtract  the  constant  contribution  of  the  numerator  of  the  first  integral 
which  removes  the  singularity  In  this  way  we  obtain: 
for  t  £  I.  2: 


Method  of  approach 
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lfv<s^s 


■VA  A  ■■■  - 


m 
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i 


hSh 

.5" 

i 

i 


i  =  J'  {[(t  -  uni  -  t  *  u)]‘  ‘fa  -  U)  -  [t(i  - 1)]1  "f(.x>] 


f:[:a  -:i]‘ 

r  t_.  fu>d; 


for  t  >  L.  2: 


dm  -;)] 

i  =  j  ■— 3— fi;)d; 

c  x 


J  L  {[(t  -  uHL  -  t  -  u)]  ‘f(t  -  u)  -  [t(L  -  t)]'  ‘f(t)| 


Results  and  discussion. 


The  present  model  w  as  applied  to  to  the  studs  of  helium  embnttlemcnt  in  fusion 
reactors  During  the  fusion  process  hebum  is  being  formed  in  the  containing  wall 
whose  embnttlemcnt  could  lead  to  senous  structural  problems  In  our  case  studs 
we  lined  up  one  or  more  hebum  bubbles  in  front  of  the  tip  of  a  typical  microcrack 
and  studied  how  ihis  affected  the  location  and  extent  of  the  plastic  zone(s)  around 
the  tip. 


In  all  of  our  calculations  we  kept  the  bubble  radius,  bubble  pressure  and  crack 
load  constant  f  or  these  quantities  we  took  the  following  typical  values  (Haubold 
19S3  and  Jaeger  et  al  19SX) 


Results  and  discussion 


bubble  radius  <2  =  0.6S.vlO  '’m 
bubble  pressure  f*  =  5ec  1 0! °.\  m 

5  •  -  *) 

CTack  load  O;  =  3x10\Ywj  ‘ 


First  we  considered  only  one  gas  bubble  at  a  distance  d  from  the  crack  tip  and 
assumed  that  onls  one  plastic  zone  is  formed:  n  (number  of  plastic  zones)  =  1. 
To  keep  a  maximum  accuracy  for  a  minimum  number  of  terms  in  our  Fourier 
expansion  \vc  set  the  crack  length  c  equal  to  d.  This  makes  the  relative  distance 
of  the  gas  bubble  from  the  crack  tip  unity. 

Fig'  5,  6  and  ~  display  the  plastic  zones  for  a  =  5S',60C  and  65‘  for  different 
values  of  d  From  these  we  see  that  the  closer  the  gas  bubble  to  the  crack  tip  the 
more  the  plastic  zone  shrinks  and  moves  away  from  the  tip  This  agrees  with  the 
global  notion  of  embnttlement  which  can  be  understood  as  a  reduction  of  the 
plastic  zone  Fig  S  shows  the  dislocation  distribution  with  and  without  bubble  for 
the  normalized  plastic  zone  They  look  quite  similar  and  generally  the  distribution 
itself  0  not  much  influenced 

So  far  we  have  considered  onls  one  dislocation  pileup  The  case  of  two  symmet- 
ncallv  placed  distributions  is  depicted  in  fig  9  Here  several  iterations  were  per¬ 
formed  The  results  point  m  the  same  direction,  namely,  the  presence  of  the  gas 
bubble  removes  and  reduces  the  plastic  zone 

Final])  the  case  of  several  equidistant  gas  bubbles  lined  up  in  front  of  the  crack 
tip  is  shown  in  fig  10  I~hc  interaction  between  the  gas  bubbles  is  negligible  for 
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our  configuration  (Ling  194S).  While  an  increase  in  the  number  of  gas  bubbles 
intensifies  the  cmbnttlemcnt  effect,  the  impact  of  each  individual  gas  bubble 
decreases  fast  as  a  function  of  its  distance  from  the  tip. 


Hence  \vc  can  conclude  that  the  presence  of  one  or  more  gas  bubbles  in  front  of 
the  crack  tip  reduces  and  removes  the  plastic  zone.  This  results  in  a  more  brittle 
behaviour  of  the  material  around  the  crack  tip  and  can  eventually  lead  to  a  brittle 
fracture. 


Summary. 

* 


The  present  analysis  has  indicated  that  the  presence  of  one  or  more  gas  bubbles 
in  front  of  a  crack  removes  and  decreases  the  plastic  zone  around  the  crack  tip. 
This  leads  to  an  embrittlement  of  the  material  and  ultimately  brittle  fracture.  The 
distance  of  the  gas  bubble  to  the  crack  tip  is  a  major  parameter  in  this  phenomenon. 
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Appendix 


We  can 
fsix.y.i) 


where 

Hi(p)  = 


express  f[x,y.i),  h{x,y.  t),  j(x,y,  i),  f{x,y,  t),  h’(x,y,  t),  hs(x,y,t)  and 
as  a  function  of  three  other  functions  H\(p),  Hi{p,  q)  and  H$(x,y,  t): 


f(x,y.t)  =  -  H.(y) 

(77) 

h(x.y.t)  =  H2(y,y) 

(78) 

j(x,y.t)  =  H2(x.y) 

(79) 

f’(x.y.t)  =  -  H|(t) 

(80) 

h’(x.y.t)  =  II2(y.t) 

(81) 

hs(x,y,t)  =  H2(x,t) 


(82) 


fs(*.y.t)  =  I  I3(x,  y.  t) 


(83) 


~  1  1  2 


V'2 


y(g]  +  6281  ) 


1  2 


R* 


*8ip  +  ?2pSi  8jp82P; 


(84) 
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—  (?!  +  ^  2^*  ^lp  +  ^lp?2B; 

*^\  2 

1  _/  -1'2 


x[-|gl”3  :e,q  +  B2qBl 


1  2 


2qgl  -  g;qg:g!  ‘]  +  ~f(gl  1  “  +  B2B;  ')  * 

2\'  2 


L  2^1  Blpq  "*"  (£2pqBl  "*"  B2p£lq  pq— 2  Blp£2q)Bl 


3  -  5  2 


I8l 


Blqglp  (glpSi  g',pg2^2Blqgl  ^ 


p,q  =  x.y  or  t 


H3(x,  y,  t)  = 


-  2x(g3t)  +  &,>.) 


(g3  +  geXg3  +  gq)2  -  4x2 


1  2 


+ 


2x(g3.  +  &jt) 


(g3  +  Be)2((g3  +  Be)2  -  4>0 


~2~  (g3y  +  &iy)(Cg3  +  6«)2  '  4x2)  + 

£  ^  L 


+  (g3  +  Be)2(B3y  +  Bey)((B3  +  Be)2  “  4*2)  ] 


in  which: 


Bi 


=  0-“  +  t' 


2  L 

-  *2) 


4xV 


g3  =  [y2  +  (t  +  x)2]’  2 
&>  =  [>2  +  (*  -  t)2]’ 2 
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